Twisted Order Parameter applied to Dimerized Ladders 
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We apply the twisted order parameter (TOP) for dimerized quantum spin ladders to locate the 
critical phases that separate gapped phases representing quantum spin liquids of various types. Using 
the DMRG, method we find that the TOP is a good order parameter for these systems regardless 
the number of legs. As a check, we reproduce with DMRG and periodic boundary conditions the 
computations previously done with Quantum Montecarlo for one-dimensional S=l/2, S=l, S=3/2 
and S=2 Heisenberg chains with alternating bonds. 
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I. INTRODUCTION 



The physics of quantum spin ladders with dimeriza- 
tion [l[ is an emblematic example of quantum phase 
transtions with a rich structure of both gapless and 
gapped phases @,[II,51 ■ The main issue in their quantum 
phase diagrams is to determine the location of critical 
points or lines, which typically separate massive phases. 
These phases exhibit quite strongly correlated proper- 
ties like Haldanc phases ||, an example of quantum 
spin liquid of great interest in condensed matter sys- 
tems. Analytical methods have been applied to uncover 
this richphase structure, some of them non-perturbative 
0, @, H B ■ However, they are not sufficient to study 
with enough precision the whole range of coupling pa- 
rameters entering in their quantum lattice Hamiltonians 
and, consequently, they have to be complemented with 
numerical studies [ToL fill . [l3j , which sometimes mod- 
ify the conclusions achieved with analytical methods, ei- 
ther quantitative and/or qualitatively. 

Quantum spin ladder systems are not mere theoreti- 
cal constructs, but they have been experimentally syn - 
thezised in several types of materials [l5|, [HI, [H, [II , 
some of them as interesting as the superconducting 
cuprates compounds or another classes of materials which 
allow for not only antiferromagnetic couplings but also 
for ferromagnetic ones as we shall be considering in this 
work. Ladder systems are also valuable candidates to im- 
plement their physics in simulable optical lattices (20ll2l|. 
where the coupling constants of the models could be ex- 
ternally manipulated more easily than with standard ma- 
terials. Furthermore, quantum spin chains and ladders 
have played a central role for ground test numerical sim- 
ulations when the density matrix renormalization group 
(DMRG) was introduced [H, 5! 51 51, 51 • 

A diversity of variants of quantum spin liquids appear 
as gapped phases in these dimerized spin ladders. These 
types of spin orders have been initially classified resorting 
to the string order parameter (SOP) [13, 53 ■ However, 
this initial SOP parameter was only valid for systems 
with spin magnitude S = 1. Then, a generalization valid 
for both integer and half-intcrger spin systems was intro- 
duced [11 and recently, we have been able to check their 



validity in a large class of dimerized spin ladders [d| and 
we have found that the generalized SOP serves as a good 
order parameter to distinguish the many quantum phases 
appearing in those ladders when their number of legs and 
spin couplings are varied. 

When moving from one gapped phase to another, this 
process generically implies crossing a critical point or line. 
Although the SOP parameter can detect the change asso- 
ciated to this phase crossing, however it does not perform 
so well when trying to locate the position of those crit- 
ical points. To this end, another order paramter called 
twisted order parameter (TOP) has been introduced for 
spin chains [11 and some ladders [3l|, 51. The TOP 
being non-local, it can be considered as a close relative 
of the SOP. In fact, we shall see that they play comple- 
mentary roles: the SOP is better suited for classifying 
quantum spin liquids, whereas the TOP is more appro- 
priate for finding their critical points. 

The idea behind the TOP comes from the twist op- 
erator introduced in the proof of the Lieb-Schulz-Mattis 
theorem (LSM) [33l . 51 in one-dimensional quantum spin 
systems. Under certain circumstances, namely, for half- 
interger S*-spin chains, the system is proved to be gap- 
less by creating a sort of twisted excited state along the 
chain, whose energy gap with respect to the ground state 
vanishes in the thermodynamic limit as 1/N, where N 
is the length of the chain. Although this twist opera- 
tor cannot be conclusive for the case of integer S'-spin 
chains, thereby opening the door for the celebrated Hal- 
dane conjecture [5|, it turns out that it can be used to 
locate the critical points of dimerized chains of arbitrary 
spins, either integer of half-integer. The way to see this 
connection is by means of the valence bond solid (VBS) 
states used as trial wavefunctions to qualitatively repre- 
sent those gapped quantum phases. 

The existence of extensions of the LSM theorem for 
spin ladders [11 can also be used as a hint to trying to 
extend the notion of the TOP parameter to more compli- 
cated lattices beyond one-dimensional systems. For spin 
ladders with S = \ spins, this generalized LSM theo- 
rem only works for an odd number of legs. Interestingly 
enough, we have found that the generalized TOP param- 
eter works for locating the critical points, regardless of 
the number of legs. 



This paper is organized as follows: in Sect|H] we in- 
troduce a set of quantum lattice Hamiltonians for spin 
ladders which present two different patterns of dimcriza- 
tion, either staggered or columnar, as well as antiferro- 
magnetic or ferromagnetic couplings between their legs 
(chains); in Sect llHI we present our numerical results for 
the TOP parameter measured in those dimcrized spin 
ladders using DMRG calculations. With this informa- 
tion we can plot the TOP vs. the dimerization strength 
and detect critical points by the vanishing of the TOP pa- 
rameter. Upon varying the dimerization patterns and the 
ferro- or antiferro-types of rung couplings, the number of 
zeros and shape of the TOP also changes in character- 
istic forms that serve us to distinguish among the set of 
quantum ladder Hamiltonians introduced in the previous 
section. ScctHVlis devoted to conclusions. 



II. MODELS OF DIMERIZED SPIN LADDERS 

In this paper we study quantum spin ladders that can 
be viewed as a certain number of bond-alternated Heisen- 
berg chains stacked one on top of another and coupled 
between them via a coupling constant J' corresponding 
to a Heiserbcrg-likc spin interaction. The bond alterna- 
tion between neighbor spins on the same leg is such that 
every strong bond is followed by a weaker one. Every leg 
can then begin either with a strong or weak coupling, and 
that initial choice determines the dimerization pattern all 
the way through the ladder. 

Then, to characterize completely our ladders we have 
to specify some parameters: the number of legs n/, the 
ferromagnetic J' < or antifcrromagnctic J' > nature 
of the Heiscnbcrg coupling between legs and the dimer- 
ization pattern of every leg, given by the coupling con- 
stant 7. 

The Hamiltonian that gives raise to this general set of 
families of dimcrized Heiscnbcrg ladders is: 

H obc = Ji,i ^ (£) • S i+1 (£) + J' Si (£) -Si(i+1) 

i i £.i 

(1) 

where Si(£) denotes a S = 1/2 spin operator at the site 
i in the leg I , with i = 1 . . . N and £ = 1 . . . m. The 
set of couplings denoted by J^i endows the ladders with 
certain dimerization patterns to be specified below. 

We shall focus our attention on two and three leg lad- 
ders, in addtition to spin chains. These cases will suf- 
fice to clarify the behavior of the TOP parameter under 
diverse dimerization patterns. Moreover, among the dif- 
ferent possibilities of establishing a dimerization pattern 
on every leg, we will only be interested in two particular 
arrangements of bonds. We will denote the staggering 
dimerization pattern that with a J^t distribution of the 
form Jij = 1 + (— 1) ?+£ 7 while the columnar pattern has 
a distribution Jij = 1 — (— 1)*7, i.e., all the legs begin 
with a strong bond. In both cases the parameter 7 is 
constrained to be within — 1 < 7 < 1. 



2 




FIG. 1: Different dimerization patterns in two and three legs 
spin ladders: a) and c) are staggered, b) and d) are columnar. 
Thicker lines correspond to a coupling constant J;^ = 1 + 7 
and thiner ones to J it i = 1 — 7 

We will also be interested in using not only open 
boundary conditions but also periodic ones. In this case, 
we have to modify the Hamiltonian fl} to take into ac- 
count the interactions between both ends: 

h pbc = h obc + J- J N/ S N (£) ■ SrW (2) 

1 

Eq. ((T|) is not translationally invariant. On the con- 
trary, if we consider ladders with an even number (greater 
than two) of sites per leg and using periodic boundary 
conditions, the Hamiltonian ^ is invariant under trans- 
lations of the form Si (£) — » Si+2 (£) and also under rever- 
sal of the dimerization strength 7 <-> —7. 

The number of possible combinations of those ingre- 
dients described above is still high, and not all of them 
may have critical properties. This fact will depend es- 
sentially on the staggering pattern and coupling constant 
J' . In particular, for ladders with two and three legs, the 
models that may have critical phases for certain config- 
urations of the microscopic parameters are the following 
ones: i/ the staggered antiferromagnetic two leg ladder, 
ii/ the columnar ferromagnetic two leg ladder, iii/ the 
staggered antifcrromagnctic three leg ladder and iv/ the 
columnar ferromagnetic three leg ladder. These models 
are depicted in fig. [TJ 

III. RESULTS 

Much like the string order parameter (SOP), the TOP 
is a non-local parameter since it involves non-local mea- 
surements. However, unlike the SOP, the success of the 
TOP relics in the 2ir twist carried out between both ends 
of the system. This fact makes it unclear how will the 
DMRG perform when measuring it. This is in sharp con- 
trast with the SOP case, which can be acuratcly mea- 
sured using only a representative subsystem within the 
bulk of the whole system. 
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FIG. 2: DMRG computation of TOP for different one dimen- 
sional Heisenberg chains with alternated bonds. The size of 
the chains is L — 64 and periodic boundary conditions have 
been used. 



A. Spin Chains 

Due to this uncertainty and as a first check, we have re- 
produced the Quantum Montecarlo results for one dimen- 
sional dimerized chains [3(| using DMRG-adapted meth- 
ods. Fig. [5] shows the TOP computed for these chains 
using periodic boundary conditions. To compute these 
values the number of DMRG sweeps was set to three, 
the number of retained states of the density matrix was 
m = 200 and the Lanczos tolerance equal to 10~ 10 . With 
these parameters, the results obtained using DMRG and 
Quantum Montecarlo agree almost identically. In order 
to go a bit further and obtain better estimates of the 
critical points using DMRG, it is more convenient to use 
open ladders. Fig. [3] represents the TOP computed us- 
ing open boundary conditions and different sizes of the 
chains. The largest size L = 150 was computed per- 
forming four DMRG sweeps, m = 300 and the tolerance 
equal to 10~ 9 , the less demanding case L = 80 required 
two sweeps, m = 120 and the tolerance equal to 10~ 8 . 
The critical points of each chain are those where the TOP 
is zero. Looking at the different graphs two observations 
arc in order: on one hand we see that the location of 
the critical points is overestimated in the case of periodic 
boundary conditions compared to the case of open chains. 
On the other hand, the values of the critical points com- 
puted with open boundary conditions are in good agree- 
ment with the values of these points computed in the 
thermodinamic limit in other works [30f . In fact, from fig. 
[3] and the curves with L = 150 we see that for the S = 1 
chain we obtain j c = 0.254, for S = 3/2 j c = 0.426. The 
case with 5 = 2 has two critical points at 7 C = 0.184 
and 7 C = 0.522. These values can be respectively com- 
pared with the previously computed "f c = 0.25997(3), 
7c = 0.43131(7),7 C = 0.1866(7),0. 5500(1). 

Regarding the critical points 7 C = in cases with 



S = 1/2 and S = 3/2, we observe that the TOP un- 
ambiguously goes to zero using periodic boundary con- 
ditions. With open boundary conditions, however, the 
TOP is still far from zero even with L = 150, although 
the tendency is clearly to vanish as we increase the size. 
Interestingly enough, a similar slow converging behaviour 
precisely at 7 C = has been also observed in [l3| mea- 
suring the string order parameter. 

B. Spin Ladders 

In the rest of this section we will measure the TOP in 
the models described previously in Sect. |TT] in order to 
test wether or not it can also be succesfully applied to 
non strictly one dimensional cases. 

We commented in the case of the spin chains that 
the TOP using open boundary conditions converged very 
slowly to zero at 7 C = in those models that in fact poss- 
esed such a critical point. Since some of the ladders that 
we will study in this section have indeed a critical point 
close to 7 = we will use periodic boundary conditions 
to avoid the slow convergence of these points. Moreover, 
as we commented before, with periodic boundary condi- 
tions our models are invariant under sign reversal of the 
dimerization strength. This means that for a given com- 
puted critical point, say 7 C! there must exist another one 
equal to —7c- The difference in the value of these oppo- 
site points will then give us a hint about the truncation 
error of our DMRG computations. 

The natural extension of the TOP as defined in [3(3 to 
systems with many legs is as follows: 

z]:=^o|exp{^^j5?}|^> (3) 
j'=i 

Where each operator now includes the contribution 
coming from each leg of the ladder Si = J2i Si (() ■ 

In fig. H] we have computed the TOP on each one of 
the critical ladders described in the previous section. 
To run the DMRG we have set m = 340, the tolerance 
equal to 10 -9 and performed three sweeps. For each one 
of the ladders we have used an arbitrary value of the 
constant J'. The critical lines of these models have been 
obtained by other techniques in references [1(3, [H, EH]. 
Since we are using periodic boundary conditions we are 
constrained to use smaller sizes than those considerered 
in the references and hence no exact agreement is 
expected. With this consideration in mind, the location 
of the critical points in fig. [4] that is, the points where 
the TOP is zero, is in good agreement with those of the 
references, as comes out from table |U 

The sign of the TOP is related with the nature of the 
VBS phase. From its definition it can be seen that its 
sign is given by a factor (— l) fc , where k is the number of 
isolated spins at the end of the chain. More explicitly, in 
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FIG. 3: TOP computed in various Heisenberg chains with 
alternated bonds and different spin values: a) S — 1/2, b) S = 
1, c) S = 3/2, d) S — 2. Computations using open boundary 
conditions properly capture the location of the critical points 
as we increase the size of the system, however the convergence 
of the TOP towards zero at the particular point y c = in the 
graphs a) and c) is very slow and makes a difference between 
both types of boundary conditions. 



a general (m, n)-VBS state, k — m if the number of sites 
is odd and k = n if it is even. 

We see in fig. [4]that the sign of the TOP is coherent with 
the VBS phases of each ladder. That is, in the cases of the 
two and three leg ladders with columnar dimcrization the 
TOP is negative from 7 = up to the corresponding crit- 
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FIG. 4: TOP computed in ladders with columnar and stag- 
gering dimerization. The constant J' has a fixed value equal 
to J' = —5 in the two and three leg ladder with colum- 
nar dimerization and J' — 1 in the two and three leg with 
staggered dimerization. Boundary conditions have been used. 
The size of the two leg ladders is L = 2 x 60, and L — 3 x 60 
for the three leg ladders 



Ladder J' 



TOP 



Other 



3, stagg. 1.0 7c = 0.001, -0.542, 0.536 -y c = 0, ±0.527 
3, col. -5.0 7 C = 0.005, -0.372, 0.368 7 C = 0, ±0.335 
2, stagg. 1.0 7c = -0.38, 0.38 7 C = ±0.34 

2, col. -5.0 7c = -0.212, 0.211 7 C = ±0.196 

TABLE I: Location of some critical points of different dimer- 
ized ladders. The identification of each ladder is done by 
means of the number of legs: two or three, and the dimeriza- 
tion pattern: staggered or columnar. We provide the values 
obtained with the TOP and the corresponding ones computed 
using different techniques. References related to these last 
values are given in the text. 



ical values 7 C where the TOP is zero, and then the TOP 
is positive all the way to 7 = 1. On the other hand the 
behavior of the TOP with the two and three leg ladders 
with staggered dimerization is the contrary. Since our 
computations have been all done using an even number 
of sites per leg, the sign of the TOP agrees with the char- 
acterization of the phases of each model done by means 
of the generalized SOP in refs. fl2l.[l3l]. That is, moving 
from 7 = to 7 = 1 the two leg ladder with columnar 
dimerization moves from a (1,1)-VBS to a (2,0)-VBS 
while the three leg ladder with columnar dimcrization 
moves from a (2, 1)-VBS to a (3,0)-VBS. In these cases 
k = 1 and the sign of the TOP is negative in the first 
phases, and k = and the sign of the TOP is positive in 
the second phases. In the case of the three leg ladder with 
staggered dimerization we move from a (1,2)-VBS to a 
(2, 1)-VBS and again the sign of the TOP is correct ac- 
cording to the value of k. We have no information on the 
VBS phases of the two leg ladder with staggered dimer- 
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ization but in this particular case the only possibility to 
obtain a negative value of the TOP is a (1, 1)-VBS, while 
a positive TOP could be a (2, 0)-VBS or a (0, 2)-VBS. 

IV. CONCLUSIONS 

The results shown in the previous section reveal that 
the DMRG algorithm is well suited to compute the 
twisted order parameter even with open boundary condi- 
tions. We have applied TOP parameter to four different 
spin ladder systems and estimated the location of their 
critical points. We find agreement with other DMRG 
calculations which compute the vanishing of the gap di- 
rectly. We have checked that the TOP computation is 
not only restricted to strictly one dimensional systems 



but it also works in two- and three-leg ladders. This 
holds true despite the LSM theorem based on the twist 
operator only works for ladders with an odd number of 
legs and S = ^ spins. In this sense, the twisted order 
parameter serves as a suitable and complementary tool 
to the string order parameter in the characterization of 
quantum phases in dimerized spin systems. 
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